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by 
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POPULAR SUMMARY 

An important component of Milankovitch’s astronomical theory of climate change 
is the precession index. The precession index, along with the Earth’s tilt and orbital 
eccentricity, are believed to he the maj or controlling factors of climate change in the lastfew 
milli on years. The precession index is e sin where e is the Earth’s orbital eccentricity 
and co s measures how close the Sun is to the Earth at midsummer. When = 90° the Sun 
is close to the Earth during northern summer, and at 270° it is far from the Earth during 
northern summer. The precession index varies with time, because both the eccentricity e 
and the parameter co^ are constantly changing due to disturbances in the Earth’s orbit by 
other planets, and due to the precession of the Earth. The change is largely periodic, with a 
period of about 23,000 years. 

This 23,000 year period is in fact observed in the paleoclimate record; deep-sea 
cores of the ooze at the bottom of the ocean clearly show that ice volume during the ice ages 
fluctuates with a 23,000 year period. (Ice volume also fluctuates with other astronomical 
periods: at 100,000 years, possibly due to changes in e alone, and at 41,000 years due to 
changes in the Earth’s tilt.) 

The standard explanation for the importance of the precession index is that cool 
northern summers are required for the growth of ice. The idea is that snow lingers through 
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the summer when the Sun is far from the Earth during northern summer, and ultimately 
builds up into ice sheets in northern Canada and Siberia. In this sense the northern 
hemisphere would be said to “control” the ice ages. 

In this paper I challenge the standard point of view. The examination of the simplest 
nonlinear energy balance climate model indicates that it is a seemingly paradoxical cooling 
of the southern oceans when the Sun’s perihelion is over the southern hemisphere (and 
therefore far from the northern hemisphere during northern summer) that causes ice to 
accumulate. This “paradoxical” cooling is an inevitable consequence of the simplest 
climate models. The cold water eventually flows north, cooling the northern hemisphere to 
the point where ice can form. In arguing for this point of view I cite evidence that the 
southern oceans cool before the northern oceans. In this sense the southern hemisphere can 
be said to “control” the precession index ice sheets. I find for the simple climate model 
considered here that the change in temperature between ice ages and temperate times can 
reach as much as 0.6° C. Presumably better climate models can increase the temperature 
difference even more. For comparison, the global average difference in temperature between 
now and the last glacial maximum was about 4° C. 
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Abstract 

A simple nonlinear energy balance climate model yields a precession index-like term 
in the temperature. Despite its importance in the geologic record, the precession index e sin 
0%, where e is the Earth’s orbital eccentricity and 6% is the Sun’s perigee in the geocentric 
frame, is not present in the insolation at the top of the atmosphere. Hence there is no one- 
for-one mapping of 23,000 and 19,000 year periodicities from the insolation to the 
paleoclimate record; a nonlinear climate model is needed to produce these long periods. A 
nonlinear energy balance climate model with radiative terms of form T a , where T is surface 
temperature and n > 1 , does produce e sin terms in temperature; the e sin (0$ terms are 

called Seversmith psychroternis. Without feedback mechanisms, the model achieves 
extreme values of ±0.64 K at the maximum orbital eccentricity of 0.06, cooling one 
hemisphere while simultaneously warming the other; the hemisphere over which perihelion 
( occurs is the cooler. In other words, the nonlinear energy balance model produces long- 

term cooling in the northern hemisphere when the Sun’s perihelion is near northern 
summer solstice and long-term warming in the northern hemisphere when the aphelion is 
near northern summer solstice. (This behavior is similar to the inertialess gray body which 
radiates like T 4 , but the amplitude is much lower for the energy balance model because of 
its thermal inertia.) This seemingly paradoxical behavior works against the standard 
Milankovitch model, which requires cool northern summers (Sun far from Earth in northern 
summer) to build up northern ice sheets, so that if the standard model is correct it must be 
more efficient than previously thought. Alternatively, the new mechanism could possibly be 
dominant and indicate southern hemisphere control of the northern ice sheets, wherein the 
southern oceans undergo a long-term cooling when the Sun is far from the Earth during 
northern summer. The cold water eventually flows north, cooling the northern hemisphere. 
This might explain why the northern oceans lag the southern ones when it comes to orbital 
■ forcing. 
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1. Introduction 


A previous paper [Rubincam, 2004] examined the simplest possible climate model: 
that of a gray body with a uniform surface albedo and infrared emissivity, but without 
thermal inertia. The gray body emits infrared radiation like T 4 , where T is surface 
temperature. It was found that the nonlinearity in T produced surface temperature terms, 
called Seversmith psychroterms, which have the form e sin a^, the same as Milankovitch’s 
precession index; here e is the Earth’s orbital eccentricity and 0% is the argument of perigee 
of the Sun’s orbit about the Earth (see Figure 1). It was suggested that the psychroterms 
were the cause of the ices ages which come and go with e sin a^, rather than Milankovitch’s 
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of examining more and more comphcated climate models to see how the psychroterms fare 
in them. 

Accordingly, the present paper examines the next simplest climate model up from 
the gray body: the nonlinear energy balance model (NEBM). The NEBM also radiates 
nonlinearly with T, but does have a thermal inertia-type term, which the gray body lacked. It 
is found here that the e sin (0^ psychroterms are still present, but reduced by an order-of- 
magnitude due to the presence of thermal inertia, to about ±0.64 K when e = 0.06. No ice- 
albedo feedback is assumed. Ice-albedo feedback might magnify the effect by a factor of 
several up to perhaps one or a few K, which would seem to be enough to produce an ice 
age. Paleoclimate data in support of the Seversmith psychroterms as the cause of the 
precession index ice ages is discussed blow. 

The precession index is an important component of Milankovitch's [1941] 
astronomical theory of climate change [e.g., Hays et al. 1976; Berger et al. 1984; Hirmov 
2000]. Once again, the precession index is e sin (0$, where e is the Earth’s orbital 
eccentricity and 0% is the argument of perigee of the Sun’s orbit about the Earth. The 
orbital geometry is given in Figure 1; note that a geocentric reference frame is assumed, so 
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that the Sun is treated as a satellite of the Earth, just as it was in two previous papers 
[Rubincam 1994, 2004]. In this frame the Sun is close to the northern hemisphere when 
= 90° and the Sun is at northern summer solstice. Likewise, the Sun is far from the 
northern hemisphere when = 270° and the Sun is at northern summer solstice. The 
precession index spectrum has a major peak at 23 kyr and a smaller one at 19 kyr (1 kyr = 
1000 years). 

It is an interesting fact that the equation for the insolation at the top of the 
atmosphere contains no terms which look like the precession index e sin a^, a result 
demonstrated by Humphreys [1964, pp. 85-87], and rediscovered by Rubincam [1994; 
1996] and Bmce G. Bills (private communication, 1994); see also Hinnov [2000] and 
Rubincam [2004], Ter ms which have the form e sin have been christened psychroterms 
[Rubincam, 2004], from the Greek word “psychro,” meaning “cold.” Thus the e sin 0) s 
term in the insolation, called the Humphreys psychroterm [Rubincam, 2004], was 
demonstrated by W. J. Humphreys to be zero. Therefore the insolation equation does not 
contain the 23 kyr and 19 kyr periods associated with e sin Because of its undoubted 
importance in the paleoclimate record, the existence of the precession index must be due to 
the Earth’s nonlinear response to the insolation. 

The standard explanation for the importance of the precession index e sin a) s is that 
cool northern summers = 270°, M s ~ 180°) allow snow to linger into the summer. 
Positive ice-albedo feedback feeds the growth of ice, ultimately building up into the major 
ice sheets [e.g., Milankovitch, 1941, pp. 435-436; Pisias andlmbrie, 1986]. This 
mechanism thus invokes a nonlinear response to obtain the 23 kyr and 19 kyr frequencies 
in ice volume [Rubincam, 1996; 2004]. Moreover, the standard “model” given by 
Milankovitch [1941] is questionable: he uses the dubious notion of caloric half-years (p. 
414), and then correlates the snow line with summer insolation and calls the correlation 
“proof’ (p. 436). 
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A previous paper [Rubincam, 2004] examined the simplest of all possible climate 
models: the gray body without thermal inertia, and where albedo A and infrared emissivity 
% are constants across its surface. The equation of the inertialess gray body is 


%0T 4 = (1-A)F s 


( 1 ) 


where T is absolute temperature, F s is the insolation and cr= 5.67 x 10' 8 Wm' 2 K 4 is the 
Stefan-Boltzmann constant. Because the gray body radiates like T 4 , it is a nonlinear climate 
model and does produce e sin 0) s terms in temperature T from the e sin (U s -free insolation. 
Any e sin ( 0 $ terms in temperature are called Seversmith psychroterms [Rubincam, 2004] to 


distinguish them from the Humphreys psychroterm in the insolation, which, as stated 


before, is zero. 

The present paper examines a slightly more complicated model: a nonlinear energy 
balance model, which radiates according to T,T 2 ,T 3 , and T 4 terms. Unlike the inertialess 
gray body examined in Rubincam [2004], the nonlinear energy balance model does have 
thermal inertia. 

In order to solve the NEBM equation (NEBM = nonlinear energy balance model), 
temperature is written a sT=T 0 + AT and T° is expanded to order (AT ) 2 . Like the gray 
body given by (1), the NEBM has a major Seversmith psychroterm term of the form AT SP «= 
-e sin 0% P,(sin <fj) in the surface temperature, with a small contribution from a similar P 3 (sin 
<p) term. Here P/sin (j>) is the Legendre polynomial of degree /, <f> is latitude, and the 


subscript “SP” on AT stands for “Seversmith psychroterm.” The magnitude of the 
temperature change can reach extremes about ±0.64 K at latitudes ±48.5° for the NEBM 
for the sum of the first two Seversmith psychroterms when the eccentricity e reaches its 
maximum value of 0.06 and co s = 90° or 270°. No ice-albedo feedback is assumed. 
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Because P,(sin 0) = sin 0 and P 3 (sin 0) = (5 sin 3 0 - 3 sin 0 )/2, both of which 
change sign when the equator is crossed, the NEBM indicates simultaneous wanning in one 
hemisphere and cooling in the other; the hemisphere over which perihelion occurs is the 
cooler. In other words, it yields precession index-type terms in T with a seemingly 
paradoxical sign: long-term cooling in the northern hemisphere when the Sun is near 
perihelion during northern summer ~ 90°) and a long-term warming in the north when 
the Sun is near aphelion during northern summer (a^ ~ 270°). The gray body does the 
same for the same reason [Rubincam, 2004]. 

The paradox is resolved by noting the following. Suppose perihelion occurs at <u^ = 
270°, so that perihelion is achieved over the southern hemisphere. Although the Sun is 
closer to the southern hemisphere with this geometry, the Sun does not spend much time 
there, because by Kepler’s second law the Sun moves fastest when it is closest to the Earth. 
On the other hand, the Sun lingers near aphelion, producing a long, cold winter in the south. 
The net result is that, on the average, the southern hemisphere is cooler. 

This effect thus works in opposition to the usual Milankovitch explanation of e sin 
0%, wherein cool summers (Sun far from Earth during northern summer) are required for 
the snow linger in order to build up into ice sheets in Canada and Fennoscandia. In this 
case the psychroterms warm the ground during the cool northern summers, making it harder 
for the standard model to work. 

The results of Rubincam [2004] and the present investigation suggest an alternative 
explanation for the e sin 6%, ice ages: the long-term cooling of the southern ocean when 0% 

~ 270° may be the trigger for the northern ice sheets. Perhaps sea ice around Antarctica 
causes ice-albedo feedback [Gildor and Tziperman, 2000, 2001] which cools the water even 
more; the cold water makes its way north, cooling the northern hemisphere and causing the 
ice sheets to grow. 

This could possibly explain the phase shift seen in the sea surface temperatures for 
the precession index, as discussed by Imbrie et al. [1988, pp. 144-148]. They note that the 
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northern oceans significantly lag the southern oceans with respect to the orbital forcing. 
This is what would be expected from the southern oceans cooling first, and then the 
northern as the cold waters spread north. Because the turn-over time of the oceans is a few 
thousand years, the northern hemisphere should lag the southern by this amount. This is 
what is observed (see the Discussion). There is some support for this point of view in the 
obliquity forcing as well; here also the northern oceans lag the southern in the obliquity 
forcing [Imbrie et al., 1988, pp. 147-148], again arguing for southern ocean control of the 
ice ages. Climate models more sophisticated than that examined here would be needed to 
pursue this idea further. 

2. Insolation 


The insolation expressed in terms of the orbital parameters is given by 


= F° 
r s 


Vfl/« (t + m)\ 


•i I ivew'i 

p ~ 0 < 7 =-<» 


COS 

sin 


l~m even 


i—m odd 


[(f-2 p)co s + (£-2p+q)M s + m(Q. s -rj-X )] (.2) 


[Rubincam, 1994, 2004], Here F s ° = 1371 W m' 2 when the reference distance r 0 is r 0 = a = 
1 AU [Hickey et al., 1988], where a is the semimajor axis of the Earth’s orbit. Also, e is 
orbital eccentricity, Q) s is the argument of perigee, M s is the mean anomaly of the Sun’s 
orbit about the Earth (see Figure 1), and is the position of the line of nodes of the Sun’s 
orbit measured from the vernal equinox; Q s = 0 hereafter. The angle £ is the Earth’s 
obliquity (currently about 23.44°), and 0 is latitude and X is east longitude on the Earth, 
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while 77 is Greenwich sidereal time. The P^fsin (f>) are the associated Legendre polynomials 

of degree £ and order m . The d t are spherical harmonic coefficients given in Rubincam 

[1994, Table 1]. The first few values are d 0 = 1/4, d x - 1/2, d^ = 5/16, and d A = -3/32. All 
the d ( with odd subscripts are zero except for d v The P’ ftnp (£) are the inclination functions 

from celestial mechanics; Rubincam [1994, Table 2; 2004, Table 2] gives them for degree 1, 
while Kaula [1966] and Caputo [1967] list them for degrees 2 through 4. The W^^e) are 

not the eccentricity functions from celestial mechanics; rather, they are special eccentricity 
functions associated solely with the insolation and are tabulated in Rubincam [1994]. 

(There is a typographical error in Table 3 of that paper; the entry for / - 2p = ±1, q = ±1 

should read 2e-3e 3 /2 instead of 2e - 3e 4 /2. Also, the (2' - 2s) ! factor in the equation 

following (7) in the same paper should read (2£- 2y)!.) 

The zonal insolation can be found from (2) by setting m = 0. The zonal insolation 
given by equation ( 8 ) in Rubincam [1994] contains some errors. The corrected equation for 
the zonal insolation at the top of the atmosphere is 

F S = Jffc)’ 

•{— /gCsin 0) • £l + — + 2ecos M + — e 2 cos2AfJ 
4 2 ^2 ^ 

+ |/>(sin 0 ) 

r e 2 

• |_(1 ) sin e sin (co + M ) + 2e sin e sin ( (O + 2 M ) 

*“) s s s s 

e 2 21 n 1 

+ — sin e sin (co — M ) + — e sin e sin (co + 3 M ) I 
8 v s s 8 s s J 
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+ — P (sin <t>) 
lo 2 


• [(i + y)(|si n2 e-j) - |(l-|e 2 )sin 2 ecos(2ffl 5 + 2Mp 
3 1 3 

+ 2e(— sin 2 £ — )cos M + — e sin 2 £ cos (2© +M ) 


9 5 3 1 

- — e sin 2 £ cos (2o) +3Af) + — e 2 (— sin 2 £ — )cos2 M 
4 5 s’ 2 4 2 5 


^•e 2 sin 2 £ cos (2fi^ ,+4M s )] +...} 


(3) 


to order e 2 . The P.(sin d) are the usual Legendre polynomials 


P o (sin0) = 1 , 
P,(sin0) = sin <p , 


3 , 1 

P 2 (sm <f>) = -sin (t> - - 


5 3 

P 3 (sin <p) = —sin 3 <j> - —sin 0 
2 2 


etc. Being just the zonal (longitude-independent) insolation, the diurnal, semidiurnal, etc. 
terms are not given in (3); but they are contained in (2) in the m * 0 terms. 

The important thing to notice about (3) is that it does not contain a term of the form 

F s = ... + coefficient x e sin CDg + .... 
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In other words, there is no term in the insolation which looks like the precession index. 
Hence there is no term in the insolation at the top of the atmosphere which has periods of 
23 kyr and 19 kyr. Bruce G. Bills, in an independent analysis (unpublished, 1994), arrived 
at the same conclusion. This lack of a precession index in the insolation was in fact 
demonstrated early on by Humphreys [1964, pp. 85-87] in his classic book on 
meteorology; Rubincam and Bills rediscovered his result via more sophisticated 
mathematics. 

The only thing in (3) which even remotely looks like the precession index is the le 
sin £ sin (0^ + 2 M s ) term which multiplies P,(sin <j>). However, this term has a high 
frequency due to the presence of the 2 M s . Its period is about half a year, not 19 kyr or 23 
kyr. 

This result has caused some confusion [j Berger , 1996; Rubincam, 1996; Hinnov , 
2000] because 23 kyr and 19 kyr periods are strongly indicated in the paleoclimate records 
[e.g.. Hays et al., 1976; Berger et al., 1984]. If the Milankovitch theory of climate is 
correct, how can there be an e sin signal in the paleoclimate records when it does not 
exist in the insolation? 

The answer is that if e sin ( 0 $ is important for climate, it must be due to the way the 
Earth responds to the insolation [Rubincam, 1994, 1996, 2004]. In other words, the 
Earth’s climate system does something nonlinear to the astronomical signal, thereby 
manufacturing the 23 kyr and 19 kyr periods. Rubincam [1994, p. 201] produced a model 
which in fact gave a precession index-like term in the radiation reaching the ground. That 
radiation is the insolation at the top of the atmosphere multiplied by the Earth’s coalbedo 
(coalbedo = 1 - A, where A is the albedo). Short-period terms in the albedo can multiply 
short-period terms* in the insolation, eliciting e sin (0$. But this “precession index” was 
extremely weak. On the other hand, the gray body elicits huge Seversmith psychroterms. 
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due the strong T 4 dependence of emitted radiation on temperature, and the tremendous 
temperature extremes caused by a lack of thermal inertia [Rubincam, 2004]. 

A simple nonlinear energy balance model (NEBM) is examined below; it produces e 
sin 6^ in the Earth’s temperature and has extreme values of about ± 0.64 K at latitudes 
±48.5° when the Earth’s orbital eccentricity is at its maximum value of e = 0.06 and = 
90° or 270°. These values are comparable to those induced by obliquity variations. No ice- 
albedo feedback is assumed. 


3. Nonlinear energy balance climate model 


The zonal nonlinear energy balance climate model (NEBM) is [e.g.. North et ai, 
1981, 1983]: 



2> ft 


D{ 

o(l 

d(i 


+ B„ + £s,r = (1 -A)F S 


1=1 


(4) 


where (i=sm<p and the B(T' terms for i > 2 are added nonlinear terms not given by North et 
al., and where now T is surface temperature in degrees Centigrade instead of the absolute 
temperature. Only the zonal terms will be needed, since the short-period diurnal terms will 
be largely damped out by the heat storage term which was not present in the inertia-free 
gray body [Rubincam, 2004] . The time t is measured in tropical years, so that (0$ + M s has 
a frequency of lit. 

The first term CcfT/dt in the equation above represents heat storage where C is a 
kind of heat capacity; thus this terms produces the thermal inertia. The second term 
represents zonal diffusion as produced by winds and ocean currents and has a diffusion 
constant D. The infrared radiation leaving the Earth in the above equation is given by 
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I = B 0 + B{T+ B 2 T 2 + B 3 T 3 + B 4 T a (5) 

The values for the B i can be estimated from the data of Graves et al. [1993] (see Figure 2); 
these values are B 0 = 195.0, 5, =1.4158, B 2 = 0.02289, B 3 = 0.001 148, and B 4 = 
0.00002089. These numbers do not come from a least-squares fit; rather, they were chosen 
to bisect the envelope of the data. 

The linearized version of (4) will be solved for first. Writing 

T=T 0 + AT (6) 


gives by (5) 


( 


I - H. + £ff,(A77 (7) 

1=1 


where 


H 0 = B 0 + B,T 0 + B 2 T 0 2 + B 3 T 0 3 + B 4 T 0 4 
H l =B l + 2B 2 T 0 + 3 B 3 T 0 2 + 4 B 4 T 0 3 
H 2 = B 2 + 3 B 3 T 0 + 6BJ 0 2 


H $ = + 4547^ 


t 


H< = B 4 
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and now T 0 = 14.9 °C. [North et al., 1981, Table 1], with H 0 = 226.005, H x = 3.1389, H 2 = 
0.1020, H 3 = 0.00235, and H 4 = 0.00002089. The linearized, first-order equation is then 


, d(T,+AT) 

dt 


Da -S)Nk±m 

dfl 


dfl 


+ H 0 + H X AT = (1 — A) 


+A K 


4 \a 


( 8 ) 


where by (3) 


AF S = P 0 (sin <j>) • \le cos M^] 

+ — P x (sin sine sin (at +M ) + 2e sine sin (co +2M )] 
2 s s s s 

+ P^fsin <$>)• [(l sin 2 e _I 

)(l + 2ecos M 


+ sin 2 e(cos (2dt)^ + M^) — 3 cos (20)^ + 3 M ,))!} (9) 


This equation is truncated at the first power in e and at /= 2. This is the usual form of the 
energy balance model and can be solved in the usual manner [North et al, 1981], Setting 

H 0 = (l-A)G 


where 
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allows the albedo A to be solved for, giving A = 0.344, in reasonable agreement with the 
observed value of 0.30 [Stephens et al., 1981], and leaving the linear equation 


, <?(A J) 
dt 


D( l-H*) 


2 Am 


d/l 


d/1 


+ + H,AT s (1 - A) AF S 


( 10 ) 


In this equation AT will be assumed to have the same form as the zonal AF S as given by (9), 
except that the parts which have a period of about a year are lagged by an angle yr : 


AT = — ecos (M s -y/) 

2 


+ — P x (sin 0) sine sin (co s + M s -y/) 
2 


5 D , • J 3 . 2 1 

+ — t, P,(sm (b)\ — sm e — 

16 2 2 U 2 


+ Tjg P 2 (sin (j)) ^ sin 2 e - ^ cos ( M s - yf) 


+ — T, eP,(sin (/>) sin 2 £ cos (2 co s + M s -yr) . (11) 

4 


The observed values are T, = 78.0 °C., t 2 = 234.8 °C., and yr= 31.5° [North et al., 
1981, Table 1]. The amplitude r, and lag angle yr were actually found by North et al. only 
for the P,( sin q>) [sin £ sin (g^ + M s -y/)] term; but they are assumed here to apply to all of 
the terms proportional to e cos (M s -yf) and e cos (2ffl^ + M s -yf) as well, since they have 
very nearly the same frequency, though this could be questioned because these terms apply 
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to different spherical harmonics. Also, North et al. [1981] calculated the value of r, by 
symmetrizing the Earth. This assumes that temperatures observed in the northern 
hemisphere also apply to the southern hemisphere. Because the continents occur mostly in 
the northern hemisphere and show greater temperature variations than the oceans, their table 
presumably overestimates the temperature extremes on the real Earth. Further, North et al. 
[1981, p. 105] assume f = 0 at the northern hemisphere’s winter solstice, while the present 
paper uses the angle 0 ) s + M s , which is zero at the vernal equinox, where the Sun is on the 
equator; thus the zero points differ by 90°. Hence in making comparisons with North et 
al .' s table, one must be careful to switch to the conventions used here. 

Equations (10) and (11) will be used to investigate the properties of the NEBM in 
the next three sections. The diffusion constant D and an associated mystery are considered 
first. 

4. Diffusion constant D 

The diffusion constant D can be found by equating the time-invariant parts of the 
P 2 (sin <p) coefficients in (10), so that 


D = 4 ( i -A)e-// 1 ^ = 0119 (12 ) 

6t 2 

The value derived (12) is much lower than the 0.649 found by North et al. [1981, p. 96] for 
their linear energy balance model; thus the present nonlinear model indicates that the overall 
heat transport is countered by other factors when averaging around the whole globe. This 
largely agrees with the result found in Rubincam [2004], where the equator-to-pole 
temperature difference on the real Earth is actually larger than that for the inertia-free gray 
body (1) with A = 0.3 and % = 1. 
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Clearly diffusion has importance for the real Earth: thanks to the Gulf Stream, 
London has a temperate climate, even though it is farther north than Winnipeg. So why is 
the diffusion constant (12) so small? 

The resolution to the mysteiy presumably lies in continentality and albedo. While 
the Gulf Stream may heat Europe, continental interiors get cold; also, the albedo increases 
with latitude [e.g., Stephens et al., 1981]. These last two would tend to counter the effects of 
winds and currents, which transport heat to the polar latitudes and moderate the equator-to- 
pole temperature difference. Thus the net effect is to give a small D in the simple model 
examined here. 

The larger North et al. [1981] value for D indicates that latitudinal heat flow is 
substantial. The discrepancy between the NEBM and the linear energy balance model may 
come from the linear model’s fitting an inherently curved set of data to the form B 0 + B X T, 
resulting in a large diffusion coefficient D, thus giving diffusion a spurious importance. 

Short et al. [1991] investigate a linear energy balance model which includes 
continentality. A similar study using a nonlinear model to investigate the discrepancy in the 
D values would be interesting and worth pursuing, but this will not be done here. 

5. Inconsistent values of C 

Finding C from (10) and (1 1) is more complicated than finding D. Writing the P,(sin 
<p) part of AT in (1 1) the usual complex notation 

f Pj (sin (j>) sin e 
and substituting in (10) yields 


2jtiCi x + 2 Dt x + H x f, 


4(1 — A) <2 
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where d) s + M s ~ 2 n (with time t measured in years), so that 


T\ = 


T x e 




where 


4(1 -A)Q 

4{2D + H x f +{2nC) 2 


(13) 


auu 


tan yf 


2jiC 
2 D+H x 


(14) 


Thus from (13) and (14) it can be seen that there are actually two equations involving C. 
Solving each equation for C using the observed values of H v Tj, and yr and the value for D 
found above yields inconsistent results: in one case C = 1.76 and C = 0.33 in the other. 
Fortunately the value for C is not needed in what follows below, and the observed values of 
Tj and y/ will simply be adopted in what comes next, rather than digress and follow up on 
the discrepancy; but it would be worth pursuing in the future. 

6. Seversmith psychroterms 


Proceeding to the second-order solution, the temperature can be written 


T= T 0 + AT + ST 


(15) 
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where 8T is the second-order part of the temperature, and the first-order part AT has already 
been given in (1 1). Substituting (15) into (10) and eliding the first-order solution leaves 


r jm 

dt 


D(l-fl 2 ) 


dm 

dfl 


dfl 


+ Hy ST + H 2 (AT) 2 
= (1 - A) AF S 


(16) 


as the second-order equation. Now the H 2 (AT) 2 term produces cross-products, where AT 
is given by (1 1). With the aid of 


P 0 (cos 0) P,(cos 0) = P,( cos 0) 

P,( cos 0) P 2 ( cos 0) = [2P,(cos 0) + 3P 3 (cos 0)]/5 


P,(cos 0) P 4 (cos 0) = [4P 3 (cos 0) + 5P 5 (cos 0)]/9 


H 2 (AT) 2 can be written 

3H 2 t 2 (e sins sin (o s ) [/j(sin 0) - — P 3 (sin 0)]/ 16 (17) 

2 

plus other terms. The important thing to note about (17) is the appearance of e sin 0%. 

Since terms of the sort (17) do not appear on the right side of (16), they must be cancelled 
by similar terms in ST in the other parts of ( 1 6). 
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The expression (17) above contains the long periods of e sin 0) s , so that Cd(SrydT 
changes slowly in (16) and can be neglected. To see this, assume ST can be written in the 
form e 2 * 1 *, where v= 271/23,000 is the frequency associated with the 23 kyr period. 
Differentiating ST with respect to time brings the v out in front of the term. Since v is small 
compared to seasonal frequencies, the time derivative is small and can be ignored. 

Hence the precession index part of the temperature can be solved for from 


^-P 3 (sin 0)3/16 


D(l-/1 2 )® 


dn\ 


dfl 


+ H l ST 5 


sp 


= -3 H 2 t , 2 {s sin e sin Q) s ) [/J (sin 0) - • 


producing 


ST sp = - e sin £ sin d) s 


IMill -. !{/> (sin 0) 
32D+16H, J l 1 


_ qd± sL Y 

2(12 D+HJ 3 Y } 


(18) 


Substituting numerical values in (18) yields 


ST s? = -13.7 e sin P^sin 0) 

+5.1 e sin 0) s P 3 (sin 0) °C 


(19) 


For e equal to the maximum value 0.06, ST SP achieves extreme values of ST sp = ± 0.64° C. 
when 0) s = 90°, 270° and (p = ± 48.5. At (p = + 48.5°, this becomes -10.67 e sin (Q s °C. 
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The results for the NEBM are graphed in Figure 3. Figure 4 illustrates the fact that the 
perihelion hemisphere is the cooler. 

It must be noted that (18) depends on the value for T p which is probably 
overestimated; see the remarks after (1 1). Hence the magnitude of the first term in (19) is 
may be too big. Also, including more terms in (11) will produce more terms in (18). 
Equation (18) is limited because (1 1) is limited to the few terms in the spherical harmonic 
expansion of temperature for which numerical values may be reasonably assumed. 

7. Comparison to obliquity changes 

The result (19) can be compared to the AT oVi expected from obliquity changes. 
From (1 1) the second harmonic in temperature is 


This will give changes in temperature as the obliquity changes with the 41 kyr period of 
magnitude by Taylor series AT ~ (dT/de)Ae 

AT obl = ^T 2 P 2 (sin0)Q-sin e cos ej Ae (20) 

where Ae is the change in obliquity. Because Ae~± l°over the obliquity cycle [e.g.. Hays 
et al. 1976; Berger et al. 1984; Hinnov 2000], 

A T obl - 40.2 Ac * ±0.7 K, (21) 

(. 
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At maximum e the Seversmith psychroterms give ± 0.64° C by (19) and are thus 90% the 
size of the obliquity variation. These calculations of course assume no feedback 
mechanisms, such as ice-albedo feedback. 

8. Temperature schematic 

Figure 5 illustrates the temperature change with and without the Seversmith 
psychroterms. In this figure 0 = +48.5° (about the same latitude as the U.S.-Canadian 
border), while e fixed at 0.06. As mentioned above, these values of (p and e are chosen to 
maximize the <5T SP in (19). 

Part (a) of Figure 5 shows the temperature variation as a function of time without the 
psychroterms. Here the temperature variation AT is represented by terms which have just 
two of the many frequencies which appear in the insolation: 

AT - 30 e cosAf s + 11.6 sin (a> s + M s ) °C. (22) 

This equation is based on (7), and omits the lag angle y/. The amplitude of the sin (cu^+ M s ) 
term, which is an annual term (equinox to equinox), comes from observing that for the 
whole Earth this term is 15.5 P,(sin <j>) sin (<o s + M s ) °C, again ignoring the lag angle 
[North et al., 1981, p. 100]; setting (/> = 448.5° gives the amplitude in the equation above. 
The e cos M s term, which cycles through 360° from perihelion to perihelion, is not observed 
separately because its period is nearly one year, making it difficult to distinguish from the 
annual term; instead by (1 1), its amplitude is scaled from the annual term, on the 
presumption that terms which have very nearly the same frequency should scale linearly 
with their importance in the insolation; however, see the caveats associated with (11) above. 

In Figure 5(a) the frequency difference between the sin (0% + M s ) and e cos M s is 
grossly exaggerated to make visible on the plot. The actual rate of change in is small; in 
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reality d) s ~ +62 arc s yr' ! [e.g.. White etal., 1996], Hence Figure 5 is schematic; the 
plotted amplitudes are realistic, but the frequencies are not, and only a few of the thousands 
of cycles are shown. The peaks show summer temperatures at this latitude, with summers 
being hottest when the Sun’s perigee 0% is near 90° and coldest near 270°, thanks to the 
presence of the e cos M s term. Similarly, the valleys show winter temperatures; winters are 
coldest at this latitude for 0% = 90° and warmest for <3% ® 270°. 

Note that there is a beat between the two terms, and that this beat period has the 
same period as e sin about 23 kyr. However, this frequency does not actually appear in 
the temperature equation above; in a signal there is a difference between the frequency of the 
envelope and the carrier wave. In the Milankovitch model, the peaks in Figure 5(a) are 
deemed to be important, and if only the peaks were plotted, a sinusoid with the 23 kyr 
period would result. However, this “strobing” at selected points in Figure 5(a) is to 
manipulate a curve which contains only short-period sin ((0$ + M s ) and e cos M s terms into 
giving the precession index frequency [Rubincam, 2004]. And to focus only on the peaks 
as important for northern ice sheet growth is to assume a nonlinear climate model; a linear 
climate model would contain only the original short-period annual and near-annual 
frequencies. 

Part (b) of Figure 5 shows the same terms as part (a), but with the addition of the 
Seversmith psychroterms discussed here. As has been shown in Section 3, for the NEBM 
this term is - 10.67 e sin °C. at the chosen latitude; the new term is displayed as the 
long-period sine curve which hugs zero. Note that its amplitude is about 1/3 that of the e 
cos M s term. Unlike the curve shown in (a), the new term has a true e sin (0$ period, so that 
there is a long-period heating and cooling. The new equation for AT graphed in (b) is thus 


AT ~30e cos M s + 1 1.6 sin (tt) s + M s ) - 10.67 e sin ® s °C. 


(23) 
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This graph is to be compared with (a). The differences are small, and are best seen 
by comparing how far the peaks rise above the AT = +10° C. line and how far the valleys 
fall below the —10° C. line. For the peaks, it will be noted that the new effect works in 
opposition to the e cos M s term, making the hot northern summers a little cooler and the 
cool summers a little hotter. Hence if the Milankovitch mechanism works, it must overcome 
the moderating influence of the new effect. 

On the other hand, examination of the valleys indicates that the Seversmith 
psychroterms make the cold northern winters even colder and the warm winters even 
warmer. Figure 4 summarizes the situation with respect to both the northern and southern 
hemispheres, and summer and winter. One could argue that these results support James 
Croll’s idea that it is very cold northern winters which cause enough ice to accumulate to 
last through the summer [e.g., Pisias and Imbrie, 1986], rather than the cool summers in the 
Milankovitch model. 

But rather than support Croll, instead it is hypothesized below in the Discussion 
section that it is the long-term hemispheric heating and cooling which might cause the ice 
sheets to grow. While Imbrie et al. [1988] argue for northern hemisphere control of the 
northern ice sheets, the effect presented here perhaps argues for southern hemisphere 
control of the northern ice sheets. The idea is that during the times when GCfe is near 270°, 
the southern hemisphere undergoes a long-term cooling, especially the southern oceans and 
perhaps the Antarctic ice cap. The cold water slowly makes its way north, cooling the whole 
Earth and eventually producing the northern ice sheets. 

9. Ice ages and the precession index 

The precession index in temperature reveals itself as the product of short-period 
terms in these (AT) 2 models. This is shown in Figure 6 in a highly schematic diagram of 
the insolation spectrum and the temperature spectrum. The Earth’s nonlinear climate 
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system creates the precession index line from the insolation. It is chiefly the P 0 (sin cp) [e 
cos Af s ] term in F s multiplying the strong annual term /^(sin cp)[ sin e sin (co s + M s )] 
which produces P,(sin cp) [e sin ta s ], giving the long periods of 23 kyr and 19 kyr. It is 
extremely interesting to note that as far as long periods are concerned, the only cross- 
product terms of any significance are these e sin 0% terms: there are no e sin 2a ^ or e to the 
first power terms, for instance. 

The magnitude of the effect is fairly large: for the energy balance model it is ±0.64 
°C. at maximum eccentricity and without ice-albedo feedback. For comparison, the 
difference in the global temperature between the present and the last glacial maximum was 
about -4° C. [Crowley, 1983, p. 868]. The obliquity oscillation gives ±0.7 °C. without 
feedback (see (21) above), so that the Seversmitn psycnroterms can get nearly as large as 
the one induced by obliquity variations. 

The sign of the new mechanism appears paradoxical: when < 2 ^ = 90° in (19), the 
Sun is close to the Earth at northern summer solstice, producing a long-term cooling in the 
northern hemisphere with simultaneous warming in the southern [Rubincam, 2004]. 
Similarly, when ( 0 $ = 270°, the Sun is far from the northern hemisphere during northern 
summer and this effect produces a long-term warming in the north, while at the same time 
the southern hemisphere undergoes long-term cooling. This is counter to intuition, which 
says that when the Sun is close to the northern hemisphere during northern summer that 
hemisphere ought to warm, not cool. However, it is the short period e cos M s term in (1 1) 
which does the warming in accordance with intuition. The peculiar sign of the mechanism 
presented here is in fact correct and necessary to achieve energy balance in the nonlinear 
models. In fact any model which radiates in the infrared like T n , where n > 1, will give a 
paradoxical e sin term. 

The sign is counter to that of the standard Milankovitch explanation of why the 
precession index is important. The standard model relies on the short period terms: it calls 
for cool northern summers when « s is near 270° and the Sun is far from the Earth, so that 
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snow lingers through the summers, ultimately building up into an ice sheet [e.g., 
Milankovitch, 1941, pp. 435-436; Pisias and Imbrie, 1986, p. 45]. By warming the 
northern cool summers somewhat and heating a little the cool winters (see Figure 4), the 
effect found here presumably makes it harder for the standard model to operate. If the 
standard explanation is responsible for the waxing and waning of ice sheets at the 23 kyr 
and 19 kyr periods, then it must be more efficient than previously thought to overcome the 
present mechanism with its opposing sign [Rubincam, 2004]. 

In contrast to smoothing out the summer temperatures, the new effect exaggerates 
the northern winter temperatures, making the cool winters colder and the warm winters 
warmer. Scottish geologist James Croll suggested in the 1860s that colder northern winters 
were the key to ice sheet growth [e.g.. Pisias and Tmhrie , 1986, p. 45]. Cfoll’s suggestion 
was later dropped in favor of Milankovitch’ s [1941] cool summers. Hence one could use 
the new effect to argue that perhaps Croll was correct, rather than Milankovitch. However, 
Croll’s case will not be argued here. 

It may be that neither the Croll nor the Milankovitch mechanisms are responsible for 
the precession index ice growth and decay; rather, it is hypothesized that the long-term 
heating and cooling of the southern hemisphere controls the northern ice sheets [Rubincam, 
2004]. This is illustrated in Figure 7. 

Figure 7 is the same as Figure 5, except that it shows temperature in the southern 
hemisphere at 0 = -48.5°. In this case the southern hemisphere undergoes a long-term 
cooling when the Sun is near <5^ = 270°, as shown in Figure 7(b) by the small but Jong 
wavelength sinusoid. The idea is that the southern oceans cool when the Sun is far from the 
Earth during northern summer. Thus, the southern oceans stay cool for a long time, 
perhaps allowing sea ice to be accumulated. The high albedo of the ice cools the southern 
hemisphere even more, so that ice-albedo feedback takes place [e.g., Gildor and Tziperman, 
2000, 2001], and the southern waters get very cold. The cold water slowly flows north. 
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eventually cooling the northern hemisphere enough to build up ice there. This hypothesis 
of course goes far beyond the simple NEBM discussed in the present paper. 

Is there observational support for the southern ocean hypothesis in the paleoclimate 
data? Imbrie et al. [1988, p. 145] show the phase lag for cold-season sea surface 
temperatures with respect to the orbital forcing for 16 sites in the North and South Atlantic 
and one in the Indian Ocean. Figure 8 (which is based on Imbrie et al. ’s figure) 
summarizes their findings, hi this diagram zero phase lag points to the top of the page, 
which corresponds to = 90°, so that the Sun is close to the northern hemisphere during 
northern summer, making the northern summers hot. The arrow labeled “MI” points to 
the phase of minimum ice volume associated with the precession index. Phase is measured 
clockwise from the top of the page, so that the minimum ice volume lags the orbital forcing 
by about 90°, according to Imbrie et al.’s convention; a 90° phase lag is what would be 
expected on a water-rich planet [e.g., Rubincam, 1995, p. 371]. Because the precession 
index has a predominantly 23 kyr period, the 90° phase lag means a time delay of about 6 
kyr. The shaded area shows the range of southern hemisphere phase lags as recovered 
from the 17 ocean-bottom core samples, while the striped area shows the range of northern 
hemisphere phase lags. There is a wide range of phase lags in the data, but they show a 
clear difference between the response times of the two hemispheres: the northern 
hemisphere significantly lags the south when it comes to the forcing. In terms of time, the 
delay is roughly 5 kyr. 

These data appear to support the hypothesis. When the Sun is close to the Earth 
during northern summer (6% ~ 90°), there is a long-term wanning in the southern 
hemisphere; thus the southern oceans should warm first, as observed. The warm waters 
flow north, wanning the northern hemisphere. Since the turn-over time of the oceans is a 
few thousand years, a time delay of a few kyr would be expected between hemispheres. The 
data give approximately 5 kyr, in apparent agreement with the hypothesis. 
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The above argument deals with warming. One would expect a similar scenario to 
hold for cooling and the growth of ice sheets, with the southern hemisphere undergoing a 
long-term cooling when the Sun is far from the Earth during northern summer ~ 270°), 
with the cold water spreading north in a few thousand years, cooling the northern 
hemisphere and allowing ice to build up in North America and Eurasia. 

There is also some support for this hypothesis in the analogous phase wheel for the 
41 kyr-period obliquity forcing, which is shown in Figure 9; this figure is based on Imbrie 
et al. ’s [1988, p. 146]. Here the phase is zero at maximum obliquity, when the poles see the 
most Sun. As in Figure 8, the diagram tends to show the southern hemisphere warming 
first, followed by the northern. The separation of phase lags is not as clear-cut as for the 
precession cycle, but this might be expected. Unlike the precession effect found here, the 
insolation and temperature variation during the obliquity cycle treats both hemispheres 
equally, as seen from (3) and (20); P 2 (sin <j>) is symmetrical about the equator. So the fact 
that there is a difference at all in the response times of the hemispheres during the obliquity 
cycle, with the southern hemisphere leading, again argues for southern hemisphere control. 
Thus there seems to be something special about the southern hemisphere, presumably 
having to do with sea ice around Antarctica; so a sea ice model [Gildor and Tziperman, 
2000; 2001] may be of some interest. 

Imbrie et al. [1988] favor northern hemisphere control of ice sheets a la 
Milankovitch. They realize that this causes a problem with regard to the precession index, 
but reject southern hemisphere forcing as “undesirably complex” [p. 154]. However, as 
shown in Figures 8 and 9, which are based on Imbrie et al .' s own diagrams, both the 
precession index and the obliquity cycle may be controlled by the southern hemisphere. 

It should also be noted that the above development was with regard to temperature 
only; no ice was included. Hence the psychroterms will exist even on an ice-free world 
[j Rubincam , 2004]. 
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10. Discussion 

Global climate models (GCMs) seek to include the greatest possible details in 
simulating climate. This comes at the expense of achieving easy physical insight, due to 
including the wealth of effects, and being forced to run numerical simulations [e.g.. North et 
al., 1981, p. 91]. In contrast, the purpose of energy balance models is to examine climate 
from the simplest possible perspective in order to gain understanding of the basics of the 
climate system. It would certainly be difficult to discover the Seversmith psychroterms 
without at least a crude analytical solution to some simple-minded model like the NEBM. 

However, energy balance models have numerous shortcomings. The heat storage 
term CoT/dt in (4) is presumably meant to be analogous to the similar term in the heat 
conduction equation pCdT/dt = KV 2 !., but it is not clear what the term means physically; 

(4) certainly does not specifically allow for heat conduction into the ground. Also, the 
relationship between surface temperature (5) and outgoing infrared radiation is empirically 
observed and not well motivated conceptually. Moreover, the diffusion constant D must be 
solved for as in (12), and by (13) and (14) the data give inconsistent values for C. 

Given these difficulties, the best approach may be to use rudimentary physical 
models with parameters whose values are already known; an example would be vertical heat 
conduction into the ground on an atmosphere-free planet with an assumed conductivity, 
density, etc. The results could then be compared to observation. Any disagreement between 
the two might indicate how to improve the model, such as to include at the next attempt 
some simple representation of an atmosphere. 

The problem with this approach is that the program may be difficult to carry out. 
Even the most elementary parameters may not be known a priori, and model complexity 
may increase rapidly. But it might be better than solving for the parameters of a poorly 
motivated model from the data: if the model is not very meaningful, what meaning can be 
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attached to the solved-for parameters? For instance, solving for D in the linear energy 
balance model gave a value that was much too large (section 4 above). 

Even assuming (4) with all its shortcomings, obstacles arose in trying to find a 
solution to the NEBM. The gray body without thermal inertia has relatively easy solution. 
It has an analytical form similar to that of the insolation (2): 






\l/4 

J 


S>,£(2 -S tm ) 


1=0 m = 0 


(l — m)! 
(£ + m)l 


•£ . M WP„(sin« 


p=0 q=-<x> 


• [(/? - 2 p)(0 s + (£-2p + q)M s + m(Q s - J] - A)] (24) 

odd 


[Rubincam, 2004], where the superscript “bb” on W^ 2p 9 (e) stands for “black body” and 
is used to distinguish these eccentricity functions from those associated with the insolation, 
which are written without the superscript. Also, the values for c t in the above equation are 

not the same as those for d t in (2). 

In contrast to the gray body, no complete analytical solution to the NEBM equation 
(4) is known. Instead, an approximation based on an iterative solution to order (AT) 2 was 
employed, and finding more than a few terms this way can be cumbersome. As a result, 
while the Seversmith psychroterms could be found to a certain level of approximation, it is 
not presently known how the P 0 (sin <j>) temperature term in the NEBM, when averaged over 
latitude, longitude, and time, behaves with e 1 . 

In the gray body (24) average temperature paradoxically goes down as e 2 increases, 
while the average insolation goes up. Specifically, if the angular brackets (7) mean 
averaging over the Earth’s surface and the bar ( T) means averaging over time, then 
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( 25 ) 


to order e 2 for the gray body, while 

W-i+f 


to the same order. The paradox is resolved by noting that 


(F s )«(r 4 )* (r) 4 , 


(26) 


(27) 


so that while the average temperature goes down as e 2 increases as shown in (25), periodic 

! 

terms are squared, cubed, and taken to the fourth power in the middle expression in (27), so 
that the power balance (26) is maintained. 

Is the same “paradoxical” behavior true of the NEBM? The answer is not known 
at present, because it is not yet clear how many iterative solutions of (4) must be found and 
how many terms must be retained before a reliable result similar to (25) can be achieved. 

Adding to the difficulties are the unknowns in (1 1); that equation was used to help 
find the psychroterms. Some of the terms in (1 1) had to be guessed at, due to a lack of 
data. 

The slow ocean currents discussed in section 9 are of course not present in the 
model discussed here, nor are continents and ice-albedo feedback [North et al, 1981, 1983; 
Graves et al . , 1993; Short et al., 1991]. The continents and feedback presumably amplify 
the effect. Ice-albedo feedback may amplify the Seversmith psychroterms to 1 or a few °C., 
which might be enough for ice sheet growth. Moreover, it may be the psychroterms which 
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triggers the ice ages which come and go with eccentricity e [Rubincam, 2004]; in other 
words, e sin (0$ — » e via some nonlinear process. These all represent avenues for future 
research. 

Because the Earth surely radiates infrared energy back into space like T a , where n > 
1 , the mechanism discussed above must certainly exist. However, the <5T SP found here may 
not be “the” precession index; some other nonlinear mechanism rather than the one 
discussed above may be stronger and be the main reason for the importance of e sin (a, for 
the ice ages. The traditional Milankovitch “model” is one possibility. The change in 
Hadley cell circulation as proposed by Lindzen [1994] and Lindzen and Pan [1994] is 
another. There are undoubtedly others. 
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Figure 1 

Geometry of the Sun’ s orbit about the Earth. Here e = obliquity, co s = argument of 
perigee, Q s = node, and f s = true anomaly. The point marked “G” on the equator marks 
the longitude of Greenwich; T| is Greenwich sidereal time measured from the vernal equinox 
VE. The subscripts “S” distinguish the Earth-based frame used here from a solar-based 
frame. Note that (0 S = 90° when perihelion occurs at northern summer solstice and (D s = 
270° when aphelion occurs at northern summer solstice. Mean anomaly M s is not shown. 
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Figure 2 

The outgoing infrared radiation as a function of temperature, based on Figure 1 of 
Graves et al. [1993]. The curve, which is a polynomial in T, is chosen to bisect the 
envelope. It is the curvature which gives rise to e sin co in temperature; a straight line would 
not do so. 
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NORTHERN HEMISPHERE 

• Cool Summers 

• Warm Winters 

• Long-term Warming 


NORTHERN HEMISPHERE 

• Hot Summers 

• Cool Winters 

• Long-term Cooling 
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Figure 4 

Summary of the long-term heating and pooling for the northern and southern 
hemispheres when cd s = 90° and 270°, due to the psychroterm effect discussed here. 
“Long-term” means on the e sin co s timescale. Note that the convention used here has CD S = 
90° when perihelion occurs at northern s umme r solstice and C0 S = 270° when aphelion 
occurs at northern summer solstice. 
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Figure 5 

Schematic temperature variation AT at (f> = 448.5°. (a) For the e cos M s and Pj(sin 
<p ) [sin (<M S 4- M s )] terms only, with e = 0.06. (b) For the e cos M s , P,( sin <j > ) [sin (fi) s 4- 
M s )], and e sin co s terms only, with e = 0.06. The e sin co s term is also shown as the small- 
amplitude sine wave. Note that the e sin term gives a long-term warming at this northern 


hemisphere location when tat - 270°. 
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Figure 6 ' 

Schematic spectrum of insolation and temperature. This diagram shows how the 
Earth’s nonlinear climate system manufactures a new spectral line in te mp erature from 
existing lines in insolation. The precession line is shown only as a single peak, and many 
other spectral lines are omitted for clarity. 
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Figure 7 

Schematic temperature variation AT at <j> = -48.5°. (a) For the e cos M s and Pj(sin 
<p ) sin (co s + M s ) terms only, with e = 0.06. (b) For the e cos M s , P,(sin <p ) [sin (co s + 
M s )l atid e sin co s terms only, with e = 0.06. The e sin 0% term is also shown as the small- 

f 

( amplitude sine wave. Note that the e sin term gives a long-term cooling at this southern 

hemisphere location when co s ~ 270°. 
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Figure 8 

Phase wheel for the precession index for cold-season sea surface temperatures and 
for ice volume for 16 sites in the North and South Atlantic and one in the Indian Ocean. 

Hot northern summers (ct) s = 90°) are at the top of the wheel at a phase lag of zero. The 
thick arrow marked “MI” is the phase lag for minimum ice volume. The phase lags for the 
southern hemisphere (SH) sites fall within the shaded area. The phase lags for the northern 
hemisphere (NH) sites fall within the striped area. The northern hemisphere clearly lags the 
southern when responding to the orbital forcing. 
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Figure 9 

Phase wheel for the obliquity cycle for cold-season sea surface temperatures and for 
ice volume for the same 17 sites as in the previous figure. Maximum obliquity is at the top 
of the wheel at a phase lag of zero. The thick arrow marked “MI” is the phase lag for 
m i n imum ice volume. The phase lags for the southern hemisphere (SH) sites fall within the 
shaded area, except for a single southern hemisphere site (marked “S” within a circle) 
which occurs near a phase lag of 190°. The phase lags for the northern hemisphere (NH) 
sites fall within the striped area. The northern hemisphere tends to lag the southern when 
responding to the orbital forcing. 
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